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We summarise applications of Dyson-Schwinger equations to the theory and 
phenomenology of hadrons. Some exact results for pseudoscalar mesons are 
highlighted with details relating to the UaW problem. We describe inferences 
from the gap equation relating to the radius of convergence for expansions of 
observables in the current-quark mass. We recapitulate upon studies of nucleon 
electromagnetic form factors, providing a comparison of the In- weighted ratios 
of Pauli and Dirac form factors for the neutron and proton. 
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1. Introduction 

The Dyson-Schwinger equations (DSEs) are a nonperturbative means of 
studying QCD in the continuum. This is illustrated in the overview pro- 
vided by Ref. [1]. A strength of the approach is that the chiral limit; viz., 
the domain of physical pion masses, is directly accessible. Indeed, with 
the established existence of a systematic, nonperturbative and symmetry- 
preserving truncation scheme [2,3] the DSEs were used to prove that the 
pion is simultaneously a Goldstone mode and a bound-state of effectively 
massive constituents [4]. Moreover, the exact nature of those dressed-quark 
and dressed-antiquark constituents is explained by the DSEs [5,6]. Namely, 
the quark-parton of QCD acquires a momentum-dependent mass function 
that at infrared momenta is larger by two orders-of-magnitude than the 
current-quark mass, an effect which owes primarily to a heavy cloud of 
gluons that clothes a low- momentum quark. a This is illustrated in Fig. 1. 



a The 2-point Schwinger functions for ghosts and gluons experience analogous nonper- 
turbative modification in the infrared [7]. 
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Fig. 1. Dressed-quark mass function, M(p): solid curves - DSE results, obtained as 
explained in Refs. [5,6], "data" — numerical simulations of unquenched lattice-QCD 
[8]. In this figure one observes the current-quark of perturbative QCD evolving into a 
constituent-quark as its momentum becomes smaller. The constituent-quark mass arises 
from a cloud of low-momentum gluons attaching themselves to the current-quark. This is 
dynamical chiral symmetry breaking: an essentially nonperturbative effect that generates 
a quark mass from nothing; viz., it occurs even in the chiral limit. 

Since a weak coupling expansion of the DSEs generates every diagram in 
perturbation theory, it is straightforward to ensure that model-dependent 
assumptions are restricted to infrared momenta. One can therefore turn a 
comparison of DSE predictions with data into a probe of the long-range 
interaction between light-quarks in QCD; namely, a unique means of ex- 
ploring light-quark confinement. In large part this study is the same as 
drawing a map of the infrared behaviour of the QCD /3-function. It is a 
fact too often ignored that the potential between infinitely heavy quarks 
measured in numerical simulations of quenched-QCD - the static potential 
- is not related in any known way to the confinement of light-quarks. 

2. Gap equation 

Since the gap equation is so fundamental to understanding hadron physics 
we reproduce it here: 



Sip)' 1 ^Z 2 ( ll -p + m hm ) + Z 1 g 2 D^(p~q) — lt ,S(q)r a „( q ,p), (1) 



where J represents a Poincare invariant regularisation of the integral, with 
A the regularisation mass-scale [4], Z? M „ is the dressed-gluon propagator, 
is the dressed-quark-gluon vertex, and m hm is the quark's A-dependent 




A' 
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bare current-mass. The vertex and quark wave-function renormalisation 
constants, Zi. 2 (( 2 , A 2 ), depend on the gauge parameter. 

The solution; namely, the dressed-quark propagator, can be written 

S(v) = 1 = Z (P 2 >t 2 ) m 

U! i~f P A(p 2 ,( 2 ) +B(p 2 ,C) ii-p + M( P 2 ) y> 

wherein the mass function, M(p 2 ), illustrated in Fig. 1, is independent of the 
renormalisation point, The solution is obtained from Eq. (1) augmented 
by the renormalisation condition S(p)~ 1 \ p2 _ ( * 2 = 17 • p + m(£ 2 ) , where 
m(C 2 ) is the running mass: Z 2 {( 2 , A 2 ) m bm (A) = Z 4 (C, 2 , A 2 ) m(( 2 ) , with 
Z± the Lagrangian-mass renormalisation constant. In QCD the chiral limit 
is strictly defined by [4]: Z 2 (( 2 ,A 2 ) m bm (A) = 0,VA 2 > ( 2 , which states 
that the renormalisation-point-invariant current-quark mass rh = 0. In this 
limit it is possible to unambiguously define the gauge invariant vacuum 
quark condensate in terms of S(p) [4,9,10]. This fact emphasises that gauge 
covariant quantities contain gauge invariant information. 

The question of whether M(p 2 ) has an expansion in current-quark mass 
around its chiral-limit value bears upon the radius of convergence for chiral 
perturbation theory. With this question one is asking whether it is possible 
to write 

oc 

M(p 2 ;m) = M(p 2 ;m = 0) + ^2a n m n (3) 



on a measurable domain of current-quark mass. It was found [10] that such 
an expansion exists for rh < m rc . The value of fh IC can be reported as fol- 
lows: for a pseudoscalar meson constituted from a quark, /, with mass rh rc 
and its antiparticle, /, m^y = 0.45 GeV. Since physical observables, such 
as the leptonic decay constant of a pseudoscalar meson, are expressed in 
terms of M(p 2 ), it follows that a chiral expansion is meaningful only for 
(mPjj,) 2 < 0.2 GeV 2 . b This entails, e.g., that it is only valid to employ chiral 
perturbation theory to fit and extrapolate results from numerical simula- 
tions of lattice-regularised QCD when the simulation parameters provide 
for m 2 < 0.2 GeV 2 . Lattice results at larger pion masses are not within the 
domain of convergence of chiral perturbation theory. 



b NB. Irrespective of the current-mass of the other constituent, a pseudoscalar meson 
containing one current-quark whose mass exceeds m cr is never within the domain of 
uniform convergence. 
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3. Mesons and the axial- vector Ward-Takahashi identity 

Dynamical chiral symmetry breaking is a fact in QCD and this amplifies 
the importance of the axial-vector Ward-Takahashi identity. The existence 
of a sensible truncation of the DSEs has enabled proof via that identity of 
a body of exact results for pseudoscalar mesons. They relate even to radial 
excitations and/or hybrids [11-13], and heavy-light [14] and heavy-heavy 
mesons [15]. The results have been illustrated using a renormalisation- 
group-improved rainbow-ladder truncation [16,17], which also provided a 
prediction of the electromagnetic pion form factor [18]. 

Implications for neutral pseudoscalar mesons have been elucidated [19]. 
In the general case the axial-vector Ward-Takahashi identity is written 

P^(k-P)=S-\k + )i 1 ,T a +i 1 ,T a S-\k_)-2iM ab T\(k-P)-A a (k-P), 

(4) 

wherein: {J ra \ a = 0, . . . , N'j — 1} are the generators of U(Nf); the dressed- 
quark propagator S = diagfS^, Sj, S s , S c , Sb, • • •] is matrix-valued; M(() is 
the matrix of rcnormalised (running) current-quark masses and Ai ab = 
tip [{F a , M}^] , where the trace is over flavour indices. The inhomoge- 
neous axial- vector vertex in Eq. (4), T^(k; P), where P is the total momen- 
tum of the quark-antiquark pair and k the relative momentum, satisfies a 
Bethc-Salpeter equation, and likewise the pseudoscalar vertex, Tl(k;P). 

The final term in Eq. (4) expresses the non-Abelian axial anomaly. It 
involves 

Au(k; P) = Jd i xd 4 ye i ^+- x - k -^N f {T° q{x) Q(0) q{yj) , (5) 

wherein the matrix element represents an operator expectation value in full 
QCD and Q(x) = e^ vpcT F^ v F^ a {x) = d^K^{x) is the topological charge 
density operator, where is the gluon field strength tensor. It is plain 
and important that only A a=0 is nonzero. 

It is fundamentally important that while Q(x) is gauge invariant, the 
associated Chern-Simons current, K^, is not. It follows that in QCD no 
physical state can couple to and hence that physical states cannot 
provide a resolution of the so-called [/^(l^problem; namely, they cannot 
play any role in ensuring that the 77' is not a Goldstone mode. 

In thinking about the C/4(l)-problem it is only necessary to focus on 
the case A a=0 ^ because if that is not so, then following Ref. [4] it is 
clear that the 7/ is certainly a Goldstone mode. A°(k; P) is a pseudoscalar 
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vertex and can therefore be expressed 



A°(k; P) = ^° 75 [%e A (k; P) + 7 ■ P^k; P) 

+ 7 • kk ■ Pg A (k; P) + a^P„H A (k; P)} . 



(6) 



Equation (4) can now be used to derive a collection of chiral-limit, pointwise 
Goldberger-Treiman relations, important amongst which is the identity 



where B (k 2 ) is obtained in solving the chiral-limit gap equation. 
It is now plain that if 



then f v 'E v i(k;0) = 0. This being true, then the homogeneous Bethe- 
Salpeter equation for the r/ docs not possess a massless solution in the 
chiral limit. The converse is also true; namely, the absence of such a solu- 
tion requires Eq. (8). Hence, Eq. (8) is a necessary and sufficient condition 
for the absence of a massless r/' bound-state. It is the chiral limit that is 
being discussed, in which case Bo(k 2 ) ^ if, and only if, chiral symme- 
try is dynamically broken. Hence, the absence of a massless r/ bound-state 
is only assured through the existence of an intimate connection between 
DCSB and an expectation value of the topological charge density. 

Reference [19] derives a range of corollaries, amongst which are mass 
formulae for neutral pseudoscalar mesons, and presents an Ansatz for the 
Bethc-Salpcter kernel that enables their illustration. Despite its simplicity, 
the model is elucidative and phenomenologically efficacious; e.g., it predicts 
77-77' mixing angles of ~ —15°, 7r°-r7 angles of ~ 1°, and a strong neutron- 
proton mass difference of 0.75 (rhd — rh u ). 

4. Baryons 

Partly because attention to this sector has only recently increased, the cur- 
rent level of expertise in studying baryons is roughly the same as it was with 
mesons more than ten years ago; viz., model building and phenomenology. 
We are a little ahead of that game, however, because much has been learnt 
in meson applications; e.g., as outlined above, a veracious understanding of 
the structure of dressed-quarks and -gluons has been acquired. 

In quantum field theory a nucleon appears as a pole in a six-point quark 
Green function. The pole's residue is proportional to the nucleon's Faddeev 
amplitude, which is obtained from a Poincare covariant Faddeev equation 



2/„,£V(fc;0) = 2B a (k 2 ) - S A (k;0) , 



(7) 



£ A (k;0) = 2B (k 2 ), 



(8) 
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that adds-up all possible quantum field theoretical exchanges and interac- 
tions that can take place between three dressed-quarks. This is important 
because modern, high-luminosity experimental facilities employ large mo- 
mentum transfer reactions; viz., Q 2 > Mfj where Mjv is the nucleon's 
mass, and hence a veracious understanding of contemporary data requires 
a Poincare covariant description of the nucleon. 

A tractable truncation of the Faddeev equation is based [20] on the ob- 
servation that an interaction which describes mesons also generates diquark 
correlations in the colour-3 channel [21]. The dominant correlations for 
ground state octet and decuplet baryons are + and 1 + diquarks because, 
e.g.: the associated mass-scales are smaller than the baryons' masses [22,23], 
namely (in GeV) - m Mo+ = 0.74 - 0.82, m (uu)i+ = m {ud)i+ = m {dd)i+ = 
0.95 — 1.02; and the electromagnetic size of these correlations is less than 
that of the proton [24] - r[ ud ] o+ « 0.7 fm, which implies r( ud ) i+ ~ 0.8 fm 
based on the p-meson/7r-meson radius-ratio [18,25]. 

The Faddeev equation's kernel is completed by specifying that the 
quarks are dressed, with two of the three dressed-quarks correlated always 
as a colour-3 diquark. Binding is then effected by the iterated exchange 
of roles between the bystander and diquark-participant quarks. A Ward- 
Takahashi-idcntity-preserving electromagnetic current for the baryon thus 
constituted is subsequently derived [26] . It depends on the electromagnetic 
properties of the axial- vector diquark correlation. 

A study of the nucleon's mass and the effect on this of a pseudoscalar 
meson cloud are detailed in [27] . Lessons learnt were employed in a series of 
studies of nucleon properties, including form factors [28-31]. The calculated 
ratio n p G p E (Q 2 )/G p M (Q 2 ) passes through zero at Q 2 w 6.5 GeV 2 [29]. The 
analogous ratio for the neutron was calculated [30]. In the neighbourhood 

2 

of Q 2 = 0, fj, n G^(Q 2 )/G r lj(Q 2 ) = — ^ Q 2 , where r n is the neutron's elec- 
tric radius. The evolution of fi p G p E (Q 2 )/G p M (Q 2 ) and ^i n G n E {Q 2 )/G n M {Q 2 ) 
on Q 2 > 2 GeV 2 are both primarily determined by the quark-core of the 
nucleon. While the proton ratio decreases uniformly on this domain [28,29], 
the neutron ratio increases steadily until Q 2 ~ 8 GeV 2 [30] . 

A form factor ratio motivated by ideas from pcrturbative QCD [32] 
is depicted in Fig. 2. The parameter A is interpreted as a mass-scale that 
defines the upper-bound on the domain of so-called soft momenta in the 
perturbative analysis. A plausible value for such a quantity is A ~ [28] . c 



C NB. A value of A ~ 0.3 GeV corresponds to a length-scale rj^ ~ lfm and it is not 
credible that perturbative QCD is applicable at ranges greater than the proton's radius. 
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Fig. 2. Weighted nucleon Pauli/Dirac form factor ratio, calculated using the framework 
of Ref. [28] and presented with A = 0.94GeV: solid curve - neutron; dashed curve - 
proton. 

It is interesting that the model of Ref. [28] yields neutron and proton ratios 
which cross at Q 2 ~ 8.5 GeV 2 . The model's prediction for truly asymptotic 
momenta is currently being explored. 

Stimulated by the possibility that Nature's fundamental "constants" 
might actually exhibit spatial and temporal variation [33], we have be- 
gun to explore the current-quark-mass-dependence of the nucleon magnetic 
moments. This complements work on hadron cr-terms [34-36]. Preliminary 
results for the quark-core contribution to this variation can be expressed 
through the following ratios evaluated at the physical current-quark mass: 



N 


p 


n 


SfiN i 5m 


0.016 


0.0042 



(9) 



It is likely that pseudoscalar meson contributions will increase these values 
by a factor of > 10 [37]. 

5. Epilogue 

Confinement and dynamical chiral symmetry breaking (DCSB) can only 
be veraciously understood in relativistic quantum field theory. The DSEs 
provide a natural vehicle for the exploration of these phenomena. 

DCSB is a singularly effective mass generating mechanism. For light- 
quarks it far outweighs the Higgs mechanism. It is understood via QCD's 
gap equation, which delivers a quark mass function with a momentum- 
dependence that connects the perturbative domain with the nonperturba- 
tive, constituent-quark domain. 
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The existence of a sensible truncation scheme enables the proof of ex- 
act results using the DSEs. That the truncation scheme is also tractable 
means the results can be illustrated, and furnishes a practical tool for the 
prediction of observables. The consequent opportunities for rapid feedback 
between experiment and theory brings within reach an intuitive under- 
standing of nonperturbative strong interaction phenomena. 

It can be argued that confinement is expressed in the analyticity prop- 
erties of dressed Schwinger functions. To build understanding it is essential 
to work toward an accurate map of the confinement force between light- 
quarks and elucidate if/how that evolves from the potential between two 
static quarks. Among the rewards are a clear connection between confine- 
ment and DCSB, an accounting of the distribution of mass within hadrons, 
and a realistic picture of hybrids and exotics. 

It is important to understand the relationship between parton properties 
on the light-front and the rest frame structure of hadrons. This is a challenge 
because, e.g., it is difficult to see how DCSB, a keystone of low-energy QCD, 
can be realised on the light- front. Parton distribution functions must be 
calculated in order to learn their content. Parametrisation is insufficient. It 
would be very interesting to know how, if at all, the distribution functions 
of a Goldstone mode differ from those of other hadrons. Answers to these 
and kindred questions are being sought using the DSEs [38,39]. 
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